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New form of matter at CERN SPS: 

Quark Matter but not Quark Gluon Plasma 

T. Csorgo 
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H - 1525 Budapest 114, POB 49 

Hungary 

I argue that a new form of matter is indeed seen in Pb+Pb collisions at CERN SPS. This Quark Matter (QM) 
is different from the theoretically predicted Quark Gluon Plasma (QGP) because its effective degrees of freedom 
seem to be the massive (dressed) constituent quarks instead of almost massless quarks and gluons. The equation 
of state of QM is hard, the time of its rehadronization is short, while the equation of state of a QGP is soft and 
the time of its rehadronization is long. Other similarities and differences are also summarized. 
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Never test for an error condition 
you don't know how to handle. " 

Steinbach's guideline 
for systems programming 



1. Introduction 

2000 has been a very exciting time in high en- 
ergy heavy ion physics, marked by two press an- 
nouncements: On 10 February, CERN summa- 
rized the results of its heavy ion program with 
claims related to the formation of a new state 
of matter in fixed target 158 AGeV Pb+Pb col- 
lisions [0J^|. On 13 July, the BNL reported on 
the observation of the first collisions of Au + Au 
nuclei at the recently completed RHIC accelera- 
tor as detected by the STAR, the PHENIX, the 
PHOBOS and BRAHMS collaborations §. 

The purpose of the heavy ion program at 
CERN SPS and at BNL RHIC accelerators is to 
produce a new state of matter, the quark gluon 
plasma (QGP), where color degrees of freedom 
are deconfined and the basic degrees of freedom 
are quarks and gluons. Due to their relatively 
large color degrees of freedom, the properties of 
QGP are dominated by that of massless gluons. 
In particular, the equation of state is soft, the en- 
tropy density is high due to the large number of 
deconfined color degrees of freedom and due to 
this reason it takes a long time for the system to 



rehadronize. Direct photons are emitted from the 
QGP by the fractionally charged quarks and the 
production of the J/tp mesons is expected to be 
suppressed due to color screening. 

Based on a critical review of the experimental 
results of the CERN heavy ion program (summa- 
rized in rcf. [0,0]), I argue that such a QGP state 
has not yet been reached in the CERN SPS heavy 
ion experiments. The circumstantial evidence 
seems to point towards the formation of a quark 
matter ( QM), where the dominant degrees of free- 
dom are massive (dressed, constituent) quarks, 
but the of effective role of gluons is secondary. 
In the followings I will argue, that in contrast 
to the soft QGP equation of state, the equation 
of state of QM has to be hard to explain the 
strong three-dimensional expansion observed in 
the hadronic final state, based on a combination 
of the data analysis in refs. f|-§]. The entropy 
density of QM seems to be relatively low ( with 
no experimental evidence for a dominant role of 
gluons). This low initial entropy density and the 
strong three-dimensional expansion of QM may 
then lead to a sudden rehadronization, where the 
abundances of directly produced particles are de- 
termined with the help of quark combinatorics as 
described in refs. |27|. Such combinatorics may 
also govern J/ip and other charmed particle pro- 
duction, a topic that deserves further investiga- 
tions 
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1.1. Summary of CERN Announcement 

On February 10, the official CERN press release 
|l],|| summarized the results of its heavy ion pro- 
gram as follows: "Compelling evidence now exists 
for the formation of a new state of matter at en- 
ergy densities at about 20 times larger than that 
in the center of atomic nuclei and temperatures 
about 100000 times higher than in the center of 
the sun. This state exhibits characteristic prop- 
erties which cannot be understood with conven- 
tional hadronic dynamics (i) , but are qualitatively 
consistent with expectations from the formation 
of a new state of matter (ii) in which quarks and 
gluons no longer feel the constraints of color con- 
finement (iii)." Both the CERN press release |l|] 
and the summary manuscript of Heinz and Ja- 
cob ^ clearly and consciously distinguished be- 
tween the claim for evidence for a new state of 
matter (that they claimed) and between the "dis- 
covery of QGP" (that they did not claim). The 
physical picture was summarized in the following 
straightforward manner: 

"(1.) Two colliding nuclei deposit energy in 
the reaction zone. The energy materializes in the 
form of quarks and gluons, which strongly inter- 
act with each other. 

(2.) This early very dense state has an en- 
ergy density of 3-4 GeV/fm 3 and the equivalent 
of a temperature of around 240 MeV. The condi- 
tions suppress the number of J-^-s (charmonia), 
enhance strangeness, and begin to drive the ex- 
pansion of the fireball. 

(3.) The "quark-gluon plasma" cools down and 
becomes more dilute. 

(4.) At an energy density of 1 GeV/fm 3 (and 
a temperature of 170-180 MeV), the quarks and 
gluons condense into hadrons, and the final abun- 
dances of the different types of particles are fixed. 

(5.) At an energy density of around 50 
MeV/fm 3 (and a temperature of 100-120 MeV) 
the hadrons stop interacting completely and the 
fireball freezes out. At this point it is expanding 
at over half the speed of light." 

It had been emphasized, that the evidence for 
the above QGP picture is circumstantial and this 
picture had been put together from many little 
observations just like a complicated jigsaw puzzle. 



2. Controversies in the QGP picture 

I think it is important to highlight some of the 
controversial points in the above QGP picture, 
that are more or less well known but may have 
not been summarized before. 

It is well known that the search for a new state 
of matter in the CERN SPS heavy ion program 
focused on two kind of QGP signatures: the early, 
penetrating probes of QGP formation (photons 
and lepton pairs that do not participate in the 
hadronic processes after their creation) and and 
the late, hadronic probes that are produced when 
the relevant kinds of strong interactions become 
negligible, due to the expansion and the related 
rarification and cooling. 

It is also well known that almost any calcula- 
tion of e.g. direct photon production or dilep- 
ton emission is sensitive to the time evolution 
and the equation of state of the system, as the 
penetrating probes are emitted from the whole 
volume integrated over the time evolution of the 
hot and dense strongly interacting matter. Hence 
the emission pattern of the penetrating probes 
depends sensitively on the time evolution of the 
temperature and density profiles, for example. 

The time evolution of the temperature or the 
local rest energy / entropy/ baryon densities, on 
the other hand, depend drastically on the dimen- 
sion of the expansion. The time evolution has to 
satisfy the boundary condition, that it ends up 
on the proper hadronic final state that has been 
determined from the analysis of the late, hadronic 
signals. The hadronic final state thus imposes a 
severe constraint on the possible time evolution 
scenarios. 

2.1. The hadronic final state 

Let us follow the strategy of backwards extrap- 
olation: It has been discovered in 1994-95 that 
the hadronic final state can be reconstructed only 
from a combined analysis of single particle spec- 
tra and two-particle Bose-Einstein correlations of 
pions and kaons |J. The hadronic final state of 
Pb + Pb collisions has been reconstructed with 
this method using the Buda-Lund hydro model in 
ref. ||. For the very different kind of experiments: 
NA44, NA49 and the preliminary WA98 data on 
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the Buda-Lund hydro model data with the Buda-Lund hydro model 



4 



central Pb+Pb collisions at CERN SPS, the aver- 
age value of the temperature as well as the value 
for of the mean transverse flow, the transverse ra- 
dius for the central and the surface temperature 
of the fireball similar fitted values were obtained, 
as well as for the mean proper-time of particle 
freeze-out and the width of the freeze-out time 
distribution . Furthermore, the mean transverse 
flow is found to be about the same as that of an 
independent analysis of ref . [0 : the mean trans- 
verse flow at the transverse radius was found to 
be (u t ) — 0.55 ± 0.06, the transverse radius pa- 
rameter was Rq = 7.1 fm, the mean freeze-out 
time after the onset of the scaling longitudinal 
expansion was found to be (tq) = 5.9 ± 0.6 fm . 
The central temperature at the mean freeze-out 
time was found to be To ~ 139 ± 6 MeV, the 
surface temperature after the particle production 
was about to end was T s ~ 85 MeV, yielding an 
average freeze-out temperature of about 100 -120 
MeV. 

The presence of a quasi-linear, Hubble -like flow 
in the final stage of Pb+Pb collisions is very well 
established, based on the observed approximately 
linear rise of the effective slope parameters of 
heavier resonances with the increasing mass: 



transverse velocity profiles 



T e ff ~ T» + m(u t y 



(1) 



The reconstructed final state of ref. |(| has 
been extrapolated backwards in time by T. S. 
Biro JtJ using exact quasi-analytic solutions of 
three-dimensional (3d) relativistic hydrodynam- 
ics assuming the presence of mixed Quark-Gluon 
Plasma - hadron gas phases, as given in Figure^. 
The result implies many things: 

i) Due to the strong 3d expansion, the system 
is able to convert a large amount of latent heat 
into hadrons in a relatively short time. Assum- 
ing a phase mixture of QGP and hadron gas, the 
QGP phase should have started to hadronize as 
early as r = 1.14 fm/c. This time scale is very 
close to the canonical guess of r» = 1 fm/c that 
is frequently used in an over-simplified Bjorken 
formula jO] to estimate the initial energy den- 
sity P,p],^2|. Bjorken's formula is based on the 
assumption of a 1 dimensional self-similar flow in 



0.6 



0.4 



<x>= -| t = 2.4 tm/c 

..1/3.; . Q 



t = 6.0 fm/c 



1 2 3 „ , 4 5 6 7 

r[fm] 



Figure 3. Biro's backwards solution of the rel- 
ativistic 3d hydrodynamics at the softest point, 
c 2 s = 0, starting from the final state give by the 
Buda-Lund fits to NA44, NA49 and WA98 data 
on single particle spectra and two-particle corre- 
lations in central Pb+Pb reactions at CERN SPS. 
From ref. M 



the beam direction: 
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which implies |T^] 

e Bj ■ = 2.8 + 0. 26GeV/fm 3 at r,- = 1.14fm/c,(3) 

in contradiction with the result of the 3d rela- 
tivistic hydrodynamical solution, where 



£3d 



= e c w 0.6GeV/fm 3 at n = 1.14fm/c (4) 



Thus the estimation of the initial energy density 
with the help of Bjorken's formula is uncertain 
and can become unreliable quantitatively (a fac- 
tor of 4 over-estimate) and qualitatively (because 
it assumes a Id instead of a 3d expansion). 

ii) Even if we assume that the initial state is 
an equilibrated, thermalized QGP consisting of 
(massless) gluons and quarks, Biro's backward ex- 
trapolation implies that the volume fraction (x) 
of QGP decreases below (x) < 1/3 within the 
first 1.3 fm/c of the 3d expansion, and the re- 
hadronization is fully completed within a time 
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period of 6 fm/c. This in turn suppresses the 
production of all the penetrating probes and im- 
plies that the signal in direct photon production, 
J I ip suppression and dilepton production must be 
very week, much weaker than signals calculated 
for an optimistic one dimensional Bjorken-type 
expansion with a long-lived rehadronization. 

Hi) Inspecting figure [3] one finds that in or- 
der to generate the strong transverse flow by the 
end of the expansion as required by the experi- 
mental data, the transverse flow has to be even 
stronger at the beginning of the rehadronization 
than at the end of this phase transition, if a soft 
equation of state is assumed. Indeed, if the pres- 
sure is constant during the rehadronization, the 
flow can only decrease, due to the radial expan- 
sion. However, the backward extrapolation de- 
pends not only on the boundary condition (the 
final state requested by the data) but also on the 
equation of state (EOS). So, one may assume that 
the hypothesis of being at the softest point is in- 
correct, Cj > 0. If the equation of state is hard, 
then the radial flow can be generated during the 
transverse expansion, perhaps even a Bjorken- 
type initial condition can be connected to the 
observed final state. The sensitivity of the ex- 
trapolation on the choice of EOS to a given final 
state has been studied recently by Schlei and col- 
laborators ||, who performed a similar but for- 
ward extrapolation using a lattice QCD inspired 
QGP equation of state and a hadron gas equation 
of state. However, the Hylander model calcula- 
tions under-estimated the transverse flow in the 
final state even when using a hard, hadron gas 
equation of state, which indicates that some of 
the observed final transverse flow may already be 
generated on the primordial level in the binary 
nucleon-nucleon collisions ||]. 

iv). Biro's full 3d hydrodynamical backwards 
solution also that the time of hand-waving argu- 
ments and order of magnitude estimates is over 
for the expansion stage of the "Little Bang" fire- 
balls. 

This observations can be summarized as fol- 
lows: a soft equation of state together with a 
zero or small primodial transverse flow is in dis- 
agreement with the boundary condition imposed 
by the hadronic final state. If one accepts the 



premiss that the strong transverse flow has not 
been present from the very beginning in the ini- 
tial state, then the strong transverse flow of the 
observed hadronic final state must have been gen- 
erated by a hard equation of state. 

2.2. A few caveats on single particle spec- 
tra and two-particle correlations 

Figure ^ indicates the M t scaling of the ef- 
fective source sizes of various particles in central 
Pb+Pb collisions at CERN SPS energies. Such a 
scaling has been predicted in an analytic calcu- 
lation of the Buda-Lund hydro model (BL-H) in 
refs. (|], with a large number of caveats that can- 
not be discussed here. See also ref. ]l3| for some 
of the details and most recent results. However, 
I would like to highlight one aspect of the BL-H 
calculation: in order to be able to calculate the 
scaling function precisely, and in order to obtain 
an approximate scaling in the transverse mass 
variable M t , the BL-H model had to assume a 
non- vanishing value for the temperature inhomo- 
geneity inside the particle emitting source. In the 
center of the plane transverse to the beam axis, 
the temperature of the source had to be slightly 
hotter than at the transverse r.m.s. radius at the 
mean freeze-out time, so that to keep the point of 
maximum emittivity close to the beam axis even 
for large values of M t . This in turn implied a sat- 
uration of the effective slope parameters T* with 
the increasing values of the particle mass, 



T* = To + m(u t )' 



To 
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R G 



1 + f (K)2 + (AT) r ) 



(5) 



(0) 



and the competition between the transverse 
flow (ut) and the transverse temperature inho- 
mogeneity {^f-) r controls the M t dependence of 
the effective source sizes as well as the flattening 
of the initial linear mass dependence of the effec- 
tive slope parameters of the single-particle spec- 
tra. Very heavy particles resolve the temperature 
inhomogeneities of the source, as the Boltzmann 
factor that characterizes their abundances focuses 
strongly their production to the hottest central 
parts, where the flow effects become limited. 
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Figure 4. M t scaling of the effective source vol- 
ume as measured by M. Murray and the NA44 
Collaboration. 



Figure 5. Slope parameters for different kind of 
particles from the NA44, NA49 and WA98 col- 
laborations. From ref. Q]. 



The central temperature is [|| To » 140 MeV, 
the flattening of the slopes sets in at about m = 
1400 MeV, which implies about 10% - 15% tem- 
perature inhomogeneity in the Pb+Pb source . 
This amount is rather small and it is in agree- 
ment with the results from a combined analy- 
sis of the single-particle spectra and the two- 
particle correlations in ref. which resulted in 
(u t ) ~ 0.55±0.06 and (^f) r ~ 0.06±0.05, which 
suggests a slope parameter of about T e ff w 315 
MeV for a particle with a mass of 1400 MeV. 

For more details on this specific point and on a 
discussion of the significance of the temperature 
inhomogeneity and the transverse flow in shap- 
ing the transverse density profile, I recommend 
ref. @. 

This 10% transverse inhomogeneity in the cen- 
tral temperature becomes rather important when 
discussing the hadron-chemical composition. In 
my opinion, this temperature inhomogeneity pro- 
hibits the association of a single temperature 
value to the hadron-chemical composition. How- 
ever, it allows for the interpretation of the ob- 
served M t scaling of HBT radii and the linearly 
rising than flattening effective slope parameters of 
the single particle spectra in a self-consistent and 



controllable manner, within the same framework. 
It is also clear, that the net baryon density is inho- 
mogeneously distributed even in the central reac- 
tion zone due to surface effects, [1b = ^B{i~x,i~y). 

Due to the importance of the inhomogeneities 
of the temperature, and the net baryon num- 
ber distributions (baryon chemical potential), I 
think that it is premature to discuss the re- 
hadronization in Pb+Pb collisions at CERN SPS 
in the framework of a complete hadrochemical 
and thermal equilibrium. The simultaneous anal- 
ysis of the two-particle correlation functions and 
the single-particle spectra indicates P,[l3| that the 
production of the observable hadrons happens in 
a relatively narrow longitudinal proper-time in- 
terval, characterized by a width of At = 1.6 ±1.5 
to be compared with the mean freeze-out proper- 
time of To = 5.9 ± 0.6 ( as measured from the 
onset of the scaling longitudinal expansion). By 
the time the particle production is over, the sur- 
face of Pb + Pb collisions cools down from 139 
MeV to about 83 MeV §,|3). It is very inter- 
esting to note, that this value is similar to the 
surface temperature of T s = 82 ± 7 MeV found in 
h+p reactions as a consequence of transverse tern- 
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perature inhomogeneities, see ref. Such 
"snow-balls" with relatively low values of surface 
temperature and a possible hotter core were re- 
ported already in S + Pb reactions in ref. || . 

Other hydro parameterizations, as reviewed in 
ref. [0, frequently neglect the effects of tem- 
perature inhomogeneities during the expansion 
and particle production stage. Energy conser- 
vation implies that the temperature cannot be 
exactly constant when particles are freezing out 
in a non-vanishing period of time from a three- 
dimensionally expanding source. Fixing the tem- 
perature to a constant in this time period, one 
finds some approximate average values of freeze- 
out temperatures in the range of Tf = 110 ± 30 
MeV. 

In the physical situation of Little Bangs, ex- 
pansion competes with the drop of the pressure 
gradients, which in turn is induced by the drop 
of the temperature on the surface. If the flow is 
small enough, a sudden drop of the temperature 
on the surface may result in a sudden drop of the 
pressure gradients on the surface, which implies 
density pile- up and a formation of a "ring of fire" , 
frequently seen in images of planetary nebulae as 
well, see ref. jl3| for further details. On the other 
hand, if the flow is strong enough, it blows away 
the material from the surface, preventing the for- 
mation of such shells of fire, and an ordinary 
expanding fireball is obtained. The former case 
seems to be realized in h + p reactions measured 
by the NA22 CERN experiment: a ring of fire 
is formed in the transverse plane due to the low 
transverse flow and due to the large temperature 
inhomogeneities. Pb+Pb collisions at CERN SPS 
are characterized with large transverse flows and 
relatively small transverse temperature inhomo- 
geneities, hence they correspond to a exploding 
fireball with a more uniform, close to Gaussian 
density distribution This result indicates 

that non-trivial time-evolution of fireball hydro- 
dynamics is intimately connected with the spatial 
inhomogeneities of the temperature and the cor- 
responding density profile. 



3. The new-old picture: Quark Matter 

If one accepts the arguments of the previous 
section about the important role of tempera- 
ture inhomogeneities in the source, the question 
arises: can we draw a phase diagram about the 
rehadronization process? Can we introduce the 
concept of a unique hadrochemical and kinetic 
temperature? In my opinion, these concepts are 
limited by the above shown variation of the local 
temperature during particle production (e.g. be- 
tween 140 to 85 MeV at the kinetic freeze-out) so 
perhaps a ball-park value can be given, but the 
non-homogeneity is important and even in ide- 
alized cases the precision of a kinetic or chemical 
freeze-out "temperature" cannot be decreased be- 
low the 20 - 30 % relative error level. Thus the 
"data points" on the beautifully drawn phase di- 
agrams e.g. in ref. |l2]| have large systematic un- 
certainties, and due to this reason I think it is 
premature to conclude at present about the sepa- 
ration of a hadrochemical freeze-out temperature 
of the order of 175 MeV from the and kinetic 
freeze-out temperature of the order of 110 MeV 
as well as about the precise value of the baryo- 
chemical potential and the temperature in various 
reactions: within 20 -30 % relative systematic er- 
ror, these values can be the same. The validity of 
the method to extract these points can be ques- 
tioned not only because it relies on the concept 
of spatially homogeneous temperature distribu- 
tions, but also because the method yields a well- 
defined value for the hadrochemical and kinetic 
freeze-out temperatures even for e + e~ reactions 
at LEP [Q. However, we know that the e + e~ 
reactions are characterized by jet production and 
non-thermal fluctuation patterns like jets within 
jets within jets etc. Hence the well defined values 
of the baryon chemical potentials and the freeze- 
out temperatures in these kind of analyzes seem 
to be more characteristic to the method than to 
the physical system under consideration. 

If the concept of the chemical freeze-out is 
maintained and embedded in a model of an ex- 
ploding fireball, a statistically acceptable x 2 fit 
to the observed hadronic abundances by Rafael- 
ski and collaborators resulted [^9| in a chemi- 
cal freeze-out temperature of the order of 140 
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MeV, within errors similar to the central kinetic 
freeze-out temperature obtained from the statis- 
tically acceptable \ 2 fits to the observed single- 
particle spectra and two-particle correlation data 
of NA44, NA49 and WA98 collaborations us- 
ing the Buda-Lund hydro parameterization ||. 
Hence a clear separation of the chemical and ther- 
mal freeze-out, (steps 4) and 5) of the Introduc- 
tion) cannot be taken for granted at present. 

4. Big Bang and Little Bang 

Table 1 briefly summarizes the similarities be- 
tween the physics of the Big Bang that resulted in 
our Universe and the physics of high energy heavy 
ion collisions or "Little Bangs" that are studied in 
the laboratory at CERN SPS and at Brookhaven 
AGS and RHIC accelerators. 

5. Strangeness enhancement 

The enhancement of strange particle produc- 
tion has been long thought to carry signals of 
QGP formation. Strangeness enhancement was 
intimately related to the question of chemical 
equilibration times in a QGP. 

In the earliest discussions of high energy heavy 
ion collisions it was assumed, that in this re- 
actions a long lived, thermally and chemically 
equilibrated QGP is formed. In such a QGP a 
large amount of ss pairs can be formed, hence 
the enhancement of strangeness production was 
proposed as a signature of QGP formation pp| . 

Later, however, it was questioned if the equilib- 
rium value of the ss number can be reached or not 
during the lifetime of a QGP phase, and it was 
shown, that the rate of the q + q — > ss reaction 
is too small to reach the equilibrium values |2l[] 
Following this observation, Rafelski and Miiller 
showed that the inclusion of the g+g — ► ss reac- 
tion a large enough strangeness production rate is 
achieved ||^] which seems to be enough to reach 
the equilibrium value of the concentration of ss 
pairs. 

The problem of the recombination of quarks 
and antiquarks into hadrons was studied by Biro 
and Zimanyi in ref . [p5[ . In this work the mecha- 
nism of hadronization was assumed to be a non- 



linear, ( quark number conserving) coalescence 
process. A few years later this model was ex- 
tended by taking into account the effects of glu- 
ons: the gluons were assumed to fission into qq 
and ss pairs, increasing the number of quarks 
and antiquarks entering into the hadronization 
process described above. Further, these au- 
thors assumed a hadrochemical evolution after 
the hadronization in order to obtain the final 
hadron numbers f24| . 

If a quark-gluon plasma state is formed in a 
heavy ion reaction, one thus expects that the 
s-quark distributions are equilibrated in a few 
fm/c, and that the production of ss quarks is 
enhanced as compared to strange particle pro- 
duction in normal hadronic interactions, because 
the following reasons: i) In a QGP, the dominant 
degrees of freedom are the gluons that can eas- 
ily enhance the strangeness content in the gluon 
fusion process g + g — > ss ; ii) at T > T c chi- 
ral symmetry is (at least partially) restored, and 
the mass of strange quarks is expected to de- 
crease to m s « 150 MeV, which implies that 
the thermal production rate is relatively high, 
N(s) cx e~ m °l T oc 1, Hi) in a baryon-rich Quark 
Gluon Plasma, the Pauli blocking of the u and 
d quarks favours (at SPS energies also) the ss 
production over the uu or dd production. Fi- 
nally, one expects that the strange quarks are 
converted into hyperons during the rehadroniza- 
tion and their abundances are thus enhanced re- 
flecting the reduced threshold of strange quark 
production in the deconfined phase. 

The question arises: how important is the dom- 
inance of gluons in the above picture? Actually, 
only the following feature matters: that the pro- 
duction threshold for a ss creation is reduced 
and that this process starts to compete with the 
production of light quark pairs. Zimanyi and 
Biro studied the problem of how quarks recom- 
bine |25[ into hadrons after the gluons have al- 
ready disappeared from the system. For clarity, 
we shall refer to this situation as to the Quark 
Matter (QM), which name does not include the 
name of gluons, in contrast to the QGP acronym. 
The kaon and hyperon enhancement has been ob- 
served to be a common feature of both QM and 
QGP |§|,||,|2|,|4|]. The question arises, what are 
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the most natural observables and can one gener- 
alize the results to non-equilibrium situations? 

I have argued that particle production in 
CERN SPS heavy ion reactions happens most 
likely in a sudden, non-equilibrium manner. 
What are the consequences of this mechanism, 
can one quantify e.g. the effect of such a mech- 
anism on the hadrochemical composition of the 
produced particles? 

Fortunately, the answer to this theoretical 
question is positive, provided that hadron pro- 
duction is characterized by the sudden recombi- 
nation of constituent quarks into hadrons, as de- 
scribed by the ALCOR model j27). Also, the ex- 
perimental results at CERN SPS Pb+Pb mea- 
surements clearly and convincingly indicated a 
huge enhancement in the production of strange 
particles Q. 

Bialas realized recently, that simple relations 
hold between various (multi-strange) antibaryon 
to baryon ratios |l5) and the Budapest group has 
proven that these relations hold not only in a lin- 
ear approximation to quark recombination, but 
also they are valid in general, if the hadron pro- 
duction happens through a sudden and complete, 
non-linear recombination of constituent quarks 
into hadrons pq] . The simplest formulation of 
such a recombinative hadron production method 
is described by the ALCOR model (ALgebraic 
COalescence for Rehadronization) p7[ . The par- 
ticle abundances are connected by the following 
simple relations |19pq]: 
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In these equations, N is the number of directly 
produced nucleons, N is the directly produced 
anti-nucleons, A is the number of directly pro- 
duced A baryons etc, so care must be taken when 
comparing with the experimental data, because 
of corrections from the resonance decays. Recent 
data from the WA97, NA44 and NA49 experi- 
ments indicate, that these relations arc satisfied 



(after resonance decays are corrected for) in cen- 
tral Pb+Pb collisions at CERN SPS. This can be 
considered as a model-independent proof |^5| , ^6| 
that constituent quark degrees of freedom are lib- 
erated in these reactions and hadron production 
proceeds via a sudden and complete coalescence of 
constituent quarks to hadrons in central Pb+Pb 
collisions at CERN SPS. 

At this point, it is natural to ask whether such 
constituent quark degrees of freedom reveal them- 
selves in p+p or p+A collisions or not? This ques- 
tion has been investigated throughoutly by the 
E910 experiment at the AGS @. The detailed 
analysis of the dependence of leading baryon pro- 
duction, strange particle production and pion 
production suggests that i) baryon stopping pro- 
ceeds through a mechanism that is different from 
energy stopping; ii) most of the energy carried 
away from precursors of energetic pions can be 
used for strange particle (e.g. A production), Hi) 
a break-up picture of the projectile to constituent 
quarks is a possible explanation of the A and K s 
production, which is very close in spirit to the 
additive quark model. 

In a QGP picture, strangeness enhancement 
happens through the large gluon density that 
create ss pairs by gluon fusion. Chiral symme- 
try restoration reduces the constituent mass of 
strange quarks from 450 MeV to about 150 MeV, 
thus reducing their production threshold . 

The question arises, do the gluonic degrees of 
freedom indeed have to play a dominant role for 
an enhancement of strange quark production? I 
think that the key requirement for strangeness 
enhancement is the reduction of the production 
threshold for ss pairs. If the quarks are con- 
fined, the even the minimal excitations require 
the presence of additional light quarks to make a 
hadron; if a Quark Matter is formed, constituent 
SS pairs can be created without the need for ad- 
ditional light quarks that also reduces the produc- 
tion threshold. Hence strangeness enhancement 
can be expected both in case of QM and in case 
of QGP formation, if the SS production treshold 
is reduced. 

A very interesting theoretical explanation of 
the process of enhanced strangeness production 
and a hard equation of state was given by Levai 
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and Heinz in ref. |3C|| . Using effective, massive 
quarks and gluons to describe the lattice QCD 
equation of state, they observed that the effec- 
tive mass of gluons is strongly increased if one ap- 
proaches from above the critical temperature. At 
the same time, the mass of quarks approached the 
constituent mass and the speed of sound remained 
rather large, c 2 s > 0.15. The increase of the effec- 
tive gluon mass with decreasing temperature may 
explain why the gluonic degrees of freedom seem 
to be less evidently required by the Pb+Pb data, 
than that of the constituent quarks |30| . 

The picture of the dominant role of con- 
stituent quarks and constituent antiquarks in the 
strangeness production, the integration of the 
gluonic degrees of freedom into an effective, con- 
fining equation of state was emphasized recently 
by the calculations of Biro, Levai and Zimanyi in 
refs. [ f31p^ ] . The necessity of a fast rehadroniza- 
tion and a three-dimensional expansion has been 
realized already in ref. J33|. 

6. Charm production and J/tp suppression 

In a QGP, color degrees of freedom are liber- 
ated, both quarks and gluons become active de- 
grees of freedom. If a cc pair is created from a 
gluonic fusion, the color interaction between these 
quarks is screened and the formation of charmo- 
nium states, in particular the creation of J/tp 
mesons, is suppressed, as predicted by Matsui and 
Satz @. 

There are a number of controversies related to 
the presentation of the NA50 data on J/tp sup- 
pression. Let me mention some of them: The 
horizontal axes of some of the figures indicates L 
a variable that is thought to be characterizing the 
length of a path that the J / tp has to travel inside 
a medium. In a three-dimensionally expending, 
rarifying and cooling fireball, such a variable is 
difficult to define, not mentioning the problem 
how to determine it experimentally. On another 
plot, the ratio measured/expected is plotted ver- 
sus the initial energy density as calculated from 
Bjorken's formula. Not only the theoretical ex- 
pectations differ from model to model, but also 
the uncertainty in the initial energy density is at 
least as big as a factor of 4. Experimental data 



points should be determined as a function of mea- 
surable quantities, and not as a function of theo- 
retical calculations. 

A particularly interesting plot has been pre- 
dicted theoretically by Kharzeev, Nardi and 
Satz |3q| : the well measurable mean transverse 
momentum (pi) of the J/tp-s in a hadron gas in- 
creases monotonically with increasing transverse 
energy, due to increased number of rescatterings 
with increasing centrality. However, if a QGP 
phase is reached in the center of the collision zone, 
the production of J/tp-s are suppressed there, 
which implies that the (pt) of the J/tp-s starts 
to decrease with increasing transverse energy E t 
produced in the experiment, after an initial rise. 
In some sense, this decrease of the (pi) of the 
J/tp-s is due to the softening of the equation of 
state when a QGP is produced. 

This plot has been compiled from the available 
NA50 data by J. Nagle in ref. Q. The result did 
not show the expected decrease of the (pi) with 
increasing E t , but followed a pattern similar to 
the prediction for a hadron gas equation of state 
|35|,[36| , see figure g. This result casts doubt on a 
J/tp suppression from a QGP in central Pb+Pb 
collisions at CERN SPS. However, the detailed 
study performed in ref. (36) indicates that the E t 
dependence of the J '/tp yield is incompatible with 
a Glauber-type calculation without initial energy 
loss, the total J/tp yield can be reproduced only 
with the help of an approximately 15 % initial 
energy loss, see fig. 0. Using this energy loss, 
the mean (p t ) transverse momenta of the J ftps 
is under-predicted. 

However, this result may be a fingerprint of the 
action of a hard equation of state, and a strong 
3-dimensional expansion in the hot hadronic mat- 
ter, that eventually enhances the (pt) of the J /ip-s 
(a transverse flow effect) in agreement with the re- 
sults of the combined analysis of the single parti- 
cle spectra and the two-particle correlation func- 
tions of pions, kaons and protons in the hadronic 
final state. 

If one assumes that quark degrees of freedom 
are liberated but gluonic degrees not, and that 
the resulting quark matter has a hard equation of 
state, the abundances of various charmed mesons 
and baryons can be calculated with the help 
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Figure 6. A systematic study of the dependence 
of the mean transverse momentum of J/tp on 
the initial energy loss as a function of the trans- 
verse energy E t of the Pb+Pb collision event . 
These NA50 datapoints (full circles with error- 
bars) are compatible with no initial energy loss. 
From rcf. pq |. The (p t ) of the J/ip-s saturates 
but does not decrease with incresing centrality 
or E t , in contrast to the prediction based on a 
Quark Gluon Plasma formation picture, given by 
Kharzeev, Nardi and Satz in ref. [BEI. 
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Figure 7. A systematic study of the dependence 
of the mean number of J/ ip mesons on the initial 
energy loss as a function of the transverse energy 
Et of the Pb+Pb collision event . These NA50 
datapoints (full circles) are compatible with a 15 
% initial energy loss. From ref. |pq|. 



of the extension of quark combinatorics to the 
charm flavor. The resulting ALCOR c model im- 
plies simple relationships between the ratios of 
( mult i) charmed antibaryon to baryon ratios, gen- 
eralizing the results presented in the section on 
strangeness 
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It is particulary interesting to note, that the 
following simple relationships are predicted by 
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and the mesonic step factors are related by a sim- 
ple relationship, 
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Quark combinatorics predicts not only the 
rates but also the slope parameters of charmed 
mesons as well, as calculated recently in refs. [j37j . 

The measurements of multi-charmed anti- 
baryon to baryon ratios , the strange/charmed 
mesonic ratios of eq. (|l^), as well as the ef- 
fective m t slopes D and J/tjj mesons can thus 
provide an important constraint and test of the 
hadronization process, and can exclude or confirm 
the possibility that charmed hadron production 
happens through quark recombination and coa- 
lescence similarly to that of the strange hadrons. 



7. Penetrating probes 



Due to lack of space, time, and expertise, I can- 
not discuss these results in detail. As I mentioned 
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in the section on the particle spectra and corre- 
lations, at CERN SPS the final state indicates a 
strong 3-dimensional expansion with a hard equa- 
tion of state, which is very efficient in reducing the 
large energy and entropy densities, which implies 
the reduction of the signal (if any) carried by the 
penetrating probes. 

7.1. Direct photons 

The experimental situation on direct photon 
production has been summarized recently in 
ref. |38| and can be briefly recapitulated as fol- 
lows: 

The relative increase of direct photons over 
photons from a conventional hadronic back- 
ground is (N^-N^ <hadr )/N^ lhadr = 12%±0.8%± 
10.9%, a non-significant value. Within the errors 
this non-significant excess is constant as a func- 
tion of multiplicity. A systematic study indicates, 
that the excess of directly produced photons in 
the most central Pb + Pb collisions (as compared 
to the hadronic event generator VENUS) is less, 
than the same excess in Pb + Nb or in Pb+Ni 
collisions [|J , see Figure ||. Indirectly, this result 
suggests the lack of significant amount of decon- 
fined, undressed light quarks in the most central 
Pb+Pb collisions at CERN SPS. 

Photon production is dominated by the 7r° — > 
7+7 decay and in the observed distributions there 
is no space for a larger than 10 % contribution of 
thermally produced photons from a Quark-Gluon 
Plasma. 

This result of the CERES/NA45 experiment 
is consistent with the picture of a strong three- 
dimensional expansion which is able to reduce the 
QGP fraction (if any) to less than 1/3 within the 
first 1.3 fm/c of rehadronization jf). 

One may expect that a long lived QGP is 
signaled with an enhanced direct photon pro- 
duction together with late freeze-out times and 
large widths of the freeze-out time distribution, 
as observed from two-particle correlation studies 
{Rout >> Rside)- Thus, the lack of direct photon 
enhancement is consistent with a non-significant 
QGP production scenario and with the experi- 
mental results on a sudden particle freeze-out, 

Pout ~ Pside- 





Figure 8. A systematic study of centrality and 
target dependence of direct photon production in 
Pb +A collisions at 158 AGeV from the WA98 
collaboration, ref. The results indicate that 
the excess of direct photons is largest in semi- 
central collisions in Pb+Pb collisions (a) and in 
collisions with smaller targets like Nb (b) and Ni 
(c) than in the most central Pb+Pb collisions. 
This result suggests that the amount of light, de- 
confined quarks is not significant in the new form 
of matter created in central Pb+Pb collisions at 
CERN SPS. 
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7.2. Dilepton production 

The CERES/NA45 experiment observed the 
inclusive e + e _ invariant mass spectra and com- 
pared the result to the expectation based on 
dilepton production from a chemically equili- 
brated, thermalized hadron gas. As compared to 
this expectation, factor of 2.6 ±0.5 ±0.6 enhance- 
ment of dilepton pairs in the 0.25 < m ee < 0.7 
interval has been reported. The enhancement was 
shown to be concentrated on the low transverse 
momentum region, 200 MeV < pf < 500 MeV, 
with almost negligible enhancement in the inter- 
val 500 MeV < p% e . The various theoretical ex- 
planations of this effect focused on hadron modi- 
fication in dense matter, as summarized recently 
in ref. ||. 

However, the results can also be interpreted so 
that the hadron gas is not in full chemical equi- 
librium. An enhancement of the low pt (rj,) rf 
and u> mesons by a factor of 5 seems to be able 
to describe the observed enhancements, compare 
the left and right panels of Fig. 6 in ref. P8[ . 
Such an enhancement has been proposed to signal 
the onset of a partial E/a(1) symmetry restora- 
tion and should be detectable with the help of 
the measurement of the intercept parameter of 
the two-pion correlation functions as a function 
of the transverse mass, mt |jjof . 

The important point of the above paragraphs 
is the following: the m t dependent production 
of u, r\ and rj' is constrained by the strength of 
the two-pion Bose-Einstein correlation function, 
and this constraint can be utilized to calibrate 
the expected number of hadronic contributions to 
dilepton decays in the range of the observed ex- 
cess. Such consistency checks between two-pion 
correlation measurements and dilepton produc- 
tion data have not yet been performed as far as I 
know. 

8. Hard Quark Matter Soft QGP 

Let me summarize the results of the CERN 
heavy ion program by slightly modifying the text 
of the original announcement (where the modi- 
fications are given in italics): "Compelling evi- 
dence now exists for the formation of a new state 
of matter at energy densities at about 5 times 



larger than that in the center of atomic nuclei... 
This state exhibits characteristic properties which 
cannot be understood with conventional hadronic 
dynamics (i) , but are qualitatively consistent with 
expectations from the formation of a new state of 
matter (ii) in which valence quarks no longer feel 
the constraints of color confinement (iii) and the 
properties of matter are dominated by that of va- 
lence quarks, the role of gluons is secondary." 

The physical picture of Quark Matter forma- 
tion can be summarized in the following straight- 
forward manner: 

(1.) Two colliding nuclei deposit energy in the 
reaction zone. The energy materializes predom- 
inantly in the form of constituent quarks, which 
strongly interact with each other. 

(2.) This early, very dense state has an energy 
density of about 1 Ge V/fm 3 and the equivalent of 
a temperature of 170 — 180 MeV. The conditions 
enhance strange quark production and the hard 
equation of state begins to drive a strongly three- 
dimensional expansion of the fireball. 

(3.) Due to the three dimensional expansion, 
the "quark matter" suddenly cools down and be- 
comes very dilute. 

(4.) The charmed, strange and light valence 
quarks recombine to form the hadrons. The al- 
lowed range for the "chemical" and the "kinetic" 
freeze-out temperatures include 140 MeV, so the 
chemical and the kinetic freeze-out happens al- 
most simultaneously. 

5.) At the time of the last interaction, the 
transverse expansion is characterized by about 
10% temperature inhomogeneities , which imply 
large inhomogeneities (edge effects) in the trans- 
verse baryon density (and the baryon chemical 
potential). At the transverse rms radius the pro- 
duced hadronic matterc expands at over half the 
speed of light." 

As mentioned, the evidence for the above 
Quark Matter picture is circumstantial and this 
picture had been put together from many little 
observations just like a complicated jigsaw puz- 
zle. This picture is significantly different from the 
also circumstantial Quark Gluon Plasma picture 
reviewed in the introduction. 
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9. Summary : Quark Matter ^ Quark 
Gluon Plasma 

A new picture is presented here, of a formation 
of a Quark Matter at CERN SPS. Quark Matter 
is a new state of matter that can be created in 
high energy heavy ion collisions. Starting from 
the reconstructed final state of Pb + Pb collisions, 
certain inconsistencies are pointed out in earlier 
attempts that tried to put together the jigsaw 
picture of these heavy ion collisions on the basis 
of the formation of a Quark Gluon Plasma. 

Constituent (valence) Quark Matter is a new 
state of matter, that is different from both or- 
dinary hadronic matter and from the much ex- 
pected Quark Gluon Plasma. 

I think that we have to keep an open eye and 
look for new, unexpected phenomena. Our com- 
munity has to be able to distinguish between the 
breakup of hadrons to a Quark Matter consist- 
ing of massive constituent quarks and between 
the complete melting of QCD matter to a Quark 
Gluon Plasma, that consists of almost massless 
quarks and gluons. 

Quark model provides a successful description 
of the bulk data on hadron spectroscopy. It is 
not surprising that the same, constituent quark 
degrees of freedom are liberated in high energy 
nuclear collisions. Note that valence gluons are 
not required in hadron spectroscopy, but valence 
(constituent) quarks are required, as they carry 
conserved quantities like charge and baryon num- 
ber. It is difficult to remove a constituent quark 
from a Quark Matter, processes like Q + Q — > Q 
are forbidden due to the conservation laws. Also, 
quarks are fermions so they cannot occupy the 
same quantum state due to Pauli blocking. I 
think these are the essential reasons why a Quark 
Matter has to be characterized by a hard equation 
of state. 

On the other hand, if a Quark Gluon Plasma 
is produced, the number of degrees of freedom 
increases drastically and the channel g + g <-> g 
opens. As the abundant gluons do not carry any 
conserved charge, their number can be changed 
relatively easily. Also, gluons are bosons so any 
number of them can occupy the same quantum 
state. I think these are the essential reasons, why 



a QGP has to be characterized by a soft equation 
of state. For a more detailed analysis of the QCD 
equation of state in terms of massive quarks and 
gluons, let me recommed the work of Levai and 
Heinz |3(J, whose phenomenological predictions 
(hard equation of state near the critical temper- 
ature etc) are in agreement with the qualitative 
analysis of the Pb+Pb data at CERN SPS. 

It seems that constituent quark degrees of free- 
dom play an important role in Pb +Pb collisions 
at CERN SPS energies, which is clearly demon- 
strated by the data on strange particle produc- 
tion. However, the liberation of gluonic degrees 
of freedom, and the softening of the equation of 
state due to the dominant role of gluons in a 
QGP is in disagreement with the available data 
at present, for example the mean transverse mo- 
mentum of the J/ips does not decrease with in- 
creasing transverse energy of the events. 

A possible interpretation of this result is that a 
Quark Matter has been created in central Pb+Pb 
collisions at CERN SPS, however, a Quark Gluon 
Plasma phase has not yet been reached. 

The discovery of QGP can be expected at 
higher initial energy densities, such as produced 
by the recently started RHIC accelerator. It is 
a very exciting challenge for RHIC to distinquish 
between various new forms of matter, for example 
between Quark Matter and Quark Gluon Plasma. 

I would like to add, that the constituent quark 
dominated Quark Matter is described in ear- 
lier publications on strangeness production under 
various, perhaps more fancy names like CQP = 
constituent quark plasma, QQP = quark anti- 
quark plasma. This picture is somewhat similar 
to the valon model as well, although the valons 
were invented to consider soft scattering prob- 
lems between quarks, as described in refs. |41{ ], 
while the constituent quarks in the present pic- 
ture emerge due to the combinatorical descrip- 
tion of hadron production which is a bound-state 
formation problem. 

In the light of the recent experimental data on 
single particle spectra, two-particle correlations, 
strange particle production, direct photon and 
dileption production as well as charmed parti- 
cle production, as summarized in refs. Jj],H the 
picture of a formation of a Quark Matter (not 
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QGP) emerges. This is in a qood agreement with 
the phenomenological analysis of the lattive QCD 
equation of state by Levai and Heinz: Instead 
of an idealized, asymptotically free quark gluon 
plasma, the lattice results for the QCD equation 
of state can be parameterized phenomenologically 
as a micture of massive quarks anti-quarks, while 
the gluons pick up a large mass near T c so that 
they stop to play a dominant role, and as a con- 
sequence the QM has hard equation of state . 
The picture of Quark Matter (not Quark Gluon 
Plasma) formation at CERN SPS can be summa- 
rized in the following manner: 

(1.) Two colliding nuclei deposit energy in the 
reaction zone. The bombarding energy breaks the 
nucleons into constituent quarks and materializes 
in the form of new constituent quark-anti-quark 
pairs. (2.) This early, very dense state has an en- 
ergy density of about 1 GeV/fm 3 and the equiv- 
alent of a temperature of more than 140 MeV. 
The conditions enhance strangeness, suppress the 
number of gluons, and the hard equation of state 
drives a strong transverse and longitudinal expan- 
sion of the fireball. 

(3.) The Quark Matter cools down fast due to 
the strong three-dimensional expansion. 

(4.) At the critical temperature of deconfine- 
mcnt, the constituent quarks suddenly recombine 
into hadrons, and the final abundances of the 
different types of directly produced particles are 
fixed. The number of charmed, strange and light 
mesons and baryons is determined from quark 
combinatorics and feed downs from the decay of 
hadronic resonances. 

(5.) The chemical and thermal freeze-out tem- 
perature distributions may be rather close to each 
other, in the range of 140 MeV. The density dis- 
tribution (baryon chemical potential) is rather 
inhomogeneous, and the central temperature at 
the mean kinetic freeze-out time decreases from 
about 140 MeV to about 85 MeV at the surface. 
At the time of the last interactions the fireball is 
expanding transversally at over half the speed of 
light." 

Table 2 summarizes briefly the main phe- 
nomenological and experimental differences be- 
tween an idealized, soft Quark Gluon Plasma, and 
the old-new, non-ideal hard Quark Matter, whose 
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Figure 9. Illustration of the difference be- 
tween Quark Matter consisting of dressed quarks, 
strange quarks and Quark Gluon Plasma which 
contains in a significant amount light and strange 
(undressed) (anti-)quarks as well as gluons . 



production in central Pb + Pb reactions at CERN 
SPS seems to be consistent with the current, ex- 
citing experimental situation. 

Figure || illustrates the phenomenological sep- 
aration of the liberation of consituent (dressed) 
quark degrees of freedom in form of Quark Mat- 
ter from the theoretically predicted liberation of 
light (undressed) quarks and gluons in very en- 
ergetic heavy ion collisions. If the gluons carry 
an effective (temperature and density dependent) 
mass, the gluonic degrees of freedom appear at 
around the T = m*(T, hb, • ••) line, which effec- 
tively separates Quark Matter from Quark Gluon 
Plasma on this illustrative example of fig. |[ Note 
however, that the "line" separating the QM from 
QGP may correspond to a crossover instead of a 
phase transition in a strict sense, similarly to that 
of the transition of a heated mono-atomic gas to a 
plasma state. At this line, the dominant degrees 
of freedom change, but the functional form of the 
equation of state may not necessarily change. 

Finally, Table 3 summarizes the 5 simple but 
different steps of a hadronizing QGP and a 
hadronizing Quark Matter. 
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Tabic 2: Properties of Quark Matter (QM) versus Quark Gluon Plasma (QGP) 



Quark Gluon Plasma 



Massless gluons and current quarks: g, q, q, s, s, 



m„ 



m T , 



m g « T c 



iq — i icq 

Id Bjorkcn expansion, softest point 
Soft equation of state 
<? s — (minimum) 
Gluons are the dominant degrees of freedom 
Long lived mixed QGP+H (Maxwell) 
T„ ~ Tk — Tp 

Rout ^ ^ Rside 
g — > SS, CC 

Strangeness suppression factor 
J/ip suppression 
(p t ) of J/tp first increases then decreases with E t 



Quark Matter 



Valence quarks: Q,Q,S,S, ... 
rriQ, m-Q, m s , m-g > T c 
3d Hubble expansion 
Hard equation of state 
Cg = 1/3 (maximum) 
Gluonic degrees of freedom frozen 
Violent explosion with sudden hadronization 
T eff =T + m(u t ) 2 

Rout — Rside 

Hadron synthesis: Quark Combinatorics 

jyn = ip/p][K/K] 3 
n ccc /n ccc = [p/ P }[D/D} 3 

(p t ) of J/ip increases with E t 



Table 3: Quark Matter (QM) versus Quark Gluon Plasma (QGP) formation: Simple steps 

Quark Matter Quark Gluon Plasma 

1 Constituent quarks are created from the bombarding energy Quarks and gluons are created from the first collisions 

2 initial energy density w 1 GeV/fm 3 w 3-4 GeV/fm 3 

3 Hard equation of state, 3d expansion Soft equation of state, Id expansion 

4 Sudden hadronization, quark combinatorics Slow hadronization, dominated by gluons 

5 Kinetic and chemical freeze-out overlaps Kinetic and chemical freeze-out well separated 



